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ABSTRACT 


A  mechanistic  model  has  been  developed  which  for  the  first  time  considers  the  effect 
of  hydrogen  entry  into  a  metal  onj^^etics  of  the  hydrogen  evolution  reaction  (h.c.r.). 
The  model  enables,  computation  me  hydrogen  surface  coverage  and  surface 
concentration;^5^^^ydrogen  adsorption,  absorption,  discharge  and  recombination  rate 
constants;  andy^iii>the  h.e.r.  coverage-dependent  transfer  coefficient,  0  and  the  exchange 
current  density  io,  from  a  knowledge  of  the  steady  state  hydrogen  permeation  current, 
cathodic  charging  current,  hydrogen  overvoltage  and  hydrogen  diffusivity.  The  model 
predicts  a  linear  relationship  between  permeation  flux  and  square-root  of  the  hydrogen 
recombination  flux  and  provides  an  analytic^  method  to  determine  the  cathodic  potential 
range  for  operation  of  a  coupled  discharge-recombination  nii^anism  of  the  h.e.r.  With 
trKxlifications  the  model  can  treat  permeation  data  for  which^  the  mechanics  of  the 
discharge  step  involve  a  (proposed)  selvedge  reaction  and^  surface  hydrogen  coverages 
are  relatively  high  as  in  the  presence  of  poisons  (e.g.,  H2S  or  As2  O3).  Some  of  the 
existing  literature  data  for  hydrogen  permeation  in  iron  and  nickel  in  acid  and  alkaline 
solutions  are  successfully  analyzed. 

Key  Words:  H  j^eimeation  model;  |urface  coverage;  h^r.  rate  constants;  ^Ivedge 
reaction;  exchange  current;  hydrogen  overvoltage;  H^abStirption  and  adsorption  rate 
constants.  I'eAr'J  ^  " 
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INTRODUCTION 

Hydrogen  is  one  of  the  most  damaging  species  in  metals,  causing  hydrogen  assisted 
cracking,  blisters,  etc.  Hydrogen  entry  into  metals  from  aqueous  solutions  has  been 
studied  extensively  by  hydrogen  permeation  experiments  of  thin  samples  using  the 
Devanathan-Stachurski  cellCl*'^.  Some  studies  carried  out  on  iron  in  acid  and  alkaline 
solutions^*^*®  reveal  a  coupled  discharge-recombination  mechanism  for  the  hydrogen 
evolution  reaction  (h.e.r.)  with  diffusion  of  absorbed  hydrogen  into  the  metal  being  rate 
controlling  for  the  permeation  process.  The  models  based  on  this  mechanism  predict  a 
square-root  relationship  between  the  hydrogen  charging  current  and  the  steady  state 
hydrogen  permeation  flux  or  currentC^X  In  all  of  these  models,  it  has  been  assumed  that  the 
hydrogen  permeation  flux  is  negligible.  It  has  also  been  assumed  that  the  hydrogen 
coverage  (0)  is  quite  low,  although  the  covwage  itself  is  unknown.  However,  especially 
when  poisons  such  as  H2S  or  AS2O3  are  present  in  the  solution,  both  the  permeation 
current  and  the  hydrogen  coverage  can  be  apprcciable(*>^).  If  the  steady  state  permeation 
current  (ico)  is  proportional  to  the  square-root  of  the  charging  current  (and  the  Tafel  slope  is 
-120  mV  decade*!  and  dT|/dlogi«o  is  -240  mV  decade*!  where  q  is  the  overpotential  for  the 
h.e.r),  the  h.ej.  is  generally  considered  to  follow  a  coupled  discharge-recombination 
mechanism.  Apart  from  these  considerations  which  have  not  produced  operative  models 
yielding  the  relevant  parameters,  e.g.,  6,  the  current  models  also  do  not  explain  the 
mecharucs  of  the  intermediate  reaction  between  adsorption  and  absorption  of  hydrogen  in 
metals,  which  is  probably  correctly  assumed  in  many  cases,  but  periiaps  not  in  all,  to  be  in 
equilibrium  at  the  cathode  surface. 

The  thrust  of  this  paper  is  to  develop  a  novel  model  for  the  electrochemical  reactions 
involving  hydrogen  evolution  and  permeaticx),  enabling  calculation  of  the  rate  constants  and 
hydrogen  coverages  for  the  first  time  using  data  obtained  from  hydrogen  permeation 
experiments.  In  this  model,  the  effect  of  hydrogen  permeation  on  the  h.e.r.  is  taken  into 


4 


t 


account  and  an  analytical  method  is  developed  to  check  the  operating  mechanism  of  the 
h.e.r.  The  model  is  tested  using  literature  data  from  hydrogen  permeation  experiments. 

THEORETICAL 

I.  Development  of  the  Basic  Model 

The  charge  transfer  (hydrogen  discharge)  reaction  can  be  represented  by 
H3O+  +  e-  +  M  =  MHads  +  H2O 

and  is  followed  by  two  reactions,  which  are  hydrogen  desorption  (or  chemical 
recombination:  2  Hads  H2)  and  hydrogen  absoiption  (permeation)  into  the  metal.  In 
Older  for  these  latter  two  reactions  to  take  place,  an  intermediate  adsorption-absorption  step 
is  involved(lO»ll).  A  quantitative  derivation  of  the  relationships  between  hydrogen  surface 
coverage,  rate  constants,  transfer  coefficient,  etc.  and  hydrogen  charging  current, 
overpotential  and  permeation  flux  is  presented  below  Avith  the  following  conditions  or 
assumptions: 

(1)  Hydrogen  discharge  involves  only  a  single  electron  transfer  reaction  and 
hydrogen  evolution  occurs  by  the  chemical  recombination  reaction  (2  Hads  = 
Hi).  It  is  also  assumed  that  t) »  RT/F  or  that  the  recombination  step  is  not  rate 
determining  so  that  the  backward  reaction  (H  oxidation)  can  be  neglected. 

Also,  Langmuir  conditions  are  assumed  for  the  hydrogen  coverage  in  the 
development  of  the  basic  model  in  this  paper,  but  the  basic  model  can  be 
modified  to  include  other  surface  coverage  conditions  as  described  elsewhere(*). 


(2)  The  intermediate  adsorption-absorption  reaction  is  in  local  equilibrium. 
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(3)  The  permeation  flux,  i«,  is  described  by  a  simple  diffusion  of  hydrogen  through 
the  metal  membrane,  with  a  finite  concentration  of  hydrogen,  Cj,  at  the 
charging  surface  and  zero  concentration  at  the  exit  (anodic)  surface. 

(4)  ic  =  ioo  +  if.  where  ic  =  charging  current  density  and  if  =  steady  state  desorption 
flux  or  hydrogen  recombination  current  density.  In  other  words,  only  steady 
state  conditions  are  considered  when  no  further  trapping  will  occur,  since  all  the 
traps  for  hydrogen  in  the  metal  will  be  saturated  by  the  end  of  the  transient  stage 
of  hydrogen  permeation.  This  assumption  of  steady  state  allows  calculation  of 
if  fiom  the  measured  ic  and  Uo  values  in  hydrogen  prameation  experiments. 

From  the  first  assumption  the  backward  reaction  of  the  h.e.r.  is  negligible  and  thus 
the  charging  cunent  is 

i,  =  Fk°aH+  C-«E*^F/RT  ^  ^  (1) 

where  k  j  =  the  rate  constant  for  the  forward  reaction;  0s  =  the  surface  coverage  by 

hydrogen;  a  =  the  transfer  coefficient;  ajj+  =  the  hydrogen  ion  activity;  ki  =  the  discharge 
rate  coefficient  =  kjajj+e'*“^®^;  =  the  equilibrium  potential  for  the  h.e.r.;  i©'  =  Fki;  i© 

=  the  exchange  current  density  =  io'  (1-0©);  0e  =  the  hydrogen  surface  coverage  at 
equilibrium;  and  a  =  F/RT  =  38.94  V-l  at  T  =  SOO^K. 

Using  Assumption  (1),  the  steady  state  desorption  flux  (hydrogen  recombination 
current)  is  given  by: 

ir  =  Fk^e]  (2) 

where  k3  =  the  recombination  rate  constant.  Using  Assumption  (3),  the  steady  state 
permeation  flux  or  current  is  given  by: 


i 


i 


ioo  *»  F(Di/L)  Cs  (3) 

where  £>i  ^  the  bulk  diffiisivity  of  hydrogen  in  the  noetal,  L  » the  thickness  of  the  metal 
membrane,  and  Cs  =  the  concentration  of  hydrogen  in  the  metal  adjacent  to  the  surface. 
From  Equation  (3), 


c.  = 


f  \ 

L 


L  “  bu 


(3a) 


where  b  =  L/(FDi). 

Considering  the  intermediate  reaction  of  hydrogen  absorption  and  adsorption,  along 
with  Assumption  (2),  the  net  rate  of  hydrogen  entry  (Uo)  is  also  given  by(^0,l  I);  - 

F(kabsds-kadsCs).  Substituting  for  Ue  from  Equation  (3)  and  rearranging,  one  obtains 


e.= 


LJ 


c,=cA" 


“•abs 


(4) 


where 

k"  = 


(4a) 


If  ^  «  kadsi  then  k"  =  =  k',  where  k'  is  usually  called  the  equilibrium  absorption- 


adsorption  constant  It  will  be  shown  later  that  the  assumption  that  the  absorption- 
adsorption  constant  is  thickness  independent  i-c.,  it  is  the  same  as  k’,  can  be  valid  for 
thick  membranes  or  low  values  of  Dj.  From  Equations  (2),  (3)  and  (4), 


6 

i.  -  k-  V4  (5) 

Equation  (5)  predicts  that  the  hydrogen  permeation  flux  (ioo)  have  a  square-root 
relationship  with  the  hydrogen  evolution  (recombination)  current,  if,  where  ir  =  ic  -  Uo  from 
Assumption  (4).  This  relationship  has  been  experimentally  observed  by  Dafft  et  al.(9).  In 
the  special  case,  when  i«, « ir  (ie.  ic  =  ir).  Equation  (5)  reduces  to being  proportional  to 
Vic.  i.e.,  the  Bockris  model. 

The  Tafel  slope  for  the  h.e.r.  and  g-^^-^-are  calculated  for  this  model.  From 
Equation  (5), 

irt  / 

L  =  KiJ'^  =  K(ic-W  (5a) 

where  K  is  a  constant. 

Taking  logarithms  of  Equation  (Sa)  and  differentiating  with  respect  to  T|, 


,  dlnd--^) 

dln(ij  .^1 _ }sL 

dq  ~  2  dri  2  dq 

Since  the  second  term  on  the  right  hand  side  of  Equation  (Sb)  is  equivalent  to 


(5b) 


^dlnij  dlnO 


20^1  ) 


,  Equation  (Sb)  becomes 


dq 


dIogC 


2-^ 

. 


dq 


dlogic 


(5c) 


Whenuo«  ic 

Soot  ~  2  Set 


(5d) 
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For  example,  when  a  =  0.5  and  9  ~  0,  Sct  =  -120  mV  decade  l  (since  from 
Equation  (1),  5^5^  =  -  - )  and  so  s»,t  =  -240  mV  decade'l,  which  is  in  exact 

agreement  with  previous  models  and  findings.  However,  if  i«.  *  ic.  Equation  (5c)  yields  Soot 

=  -120  mV  decade*l  for  Set  =  -120  mV  decade**.  Thus,  in  general,  the  following 

dri 

inequality  will  hold  for 


dq 

dq 

<2 

dq 

dlogic 

< 

dlogL 

diogic 

(5c) 


It  may  be  additionally  noted  &om  Equation  (5c)  that  for  significant  hydrogen 
penneation  (Uo),  either  the  Tafel  slope  or  diVdlogioo  or  both  will  vary  as  a  consequence  of 
the  dependence  of  0  on  Tl  implicit  in  Equation  (5).  Many  data  relating  to  hydrogen 
permeation  including  those  of  Bockiis  et  al.(?\  Kato  et  alP)  and  Dafft  et  clearly  show 
the  nonconstant  behavior  of  these  quantities.  It  will  be  cautioned  here  that  an  increase  in 
ioq^  corresponds  to  an  increase  in  6$  in  which  case  the  Frumkin-Temkin  correction  may 
need  to  be  applied  for  ic  and  if  in  Equations  (1)  and  (2).  This  will  be  discussed  towards 
the  end  of  this  analysis. 

In  order  to  find  the  values  of  k",  k3  and  i’o,  analysis  of  the  Ti-ic-i-i  data  has  to  be 
carried  out  Rearranging  Equation  (1), 


es  =  i-i^ 


(6) 


Equation  (6)  can  be  called  the  polarized  adsorption  isothenn.  This  is  to  be  contrasted  with 
the  Langmuir  adsorption  isotherm  which  is  good  only  if  at  least  quasi-equilibrium  exists  so 
that  the  forward  and  the  reverse  currents  are  equal.  On  the  other  hand.  Equation  (6)  gives 


the  adsorption  isotherm  (i.e.  surface  coverage)  for  a  non-reversible  reaction,  such  as  when 
the  nonreversible  reaction  is  coupled  to  a  subsequent  reaction.  Also  from  Equations  (3a) 
and  (4), 


0s  =  ^  (6a) 

Equation  (6a)  can  be  called  the  polarized  adsorption-absorption  isotherm.  By 
differentiating  Equation  (6a)  with  respect  to  ioo. 


(6b) 


Equati(»i  (6b)  tells  us  that  the  gradient  of  the  polarized  adsorption  isotherm  (ie.  the 
polarized  coverage)  with  the  absorption  (ot  permeation)  current  will  yield  the  thickness 
dependent  absorption-adsorption  constant,  k".  However,  in  practice,  0s  is  not  known.  So 
one  has  to  eliminate  0s  between  Equations  (6)  and  (6a),  in  order  to  obtain  the  following 
relationship: 

ic  =  io’  [l-(bijk")]e-aan  (6c) 

i.e., 

icc««l  = -(bio'/k")(u,)  +  io'  (6d) 

Equation  (6d)  predicts  that  a  linear  relationship  should  be  obtained  between  ieo  and 
fj^  =  (ice®*”!)  where  a  is  obtained  by  the  procedure  shown  in  the  Appendix,  and  the  value 
of  the  bulk  hydrogen  diffusivity,  Di,  in  the  constant,  b,  is  separately  obtained  from  other 
(transient)  permeation  measurements  or  from  the  literature.  Thus,  Equation  (6d)  enables 
direct  evaluation  of  k"  (from  the  slope  of  fit^  vs  i«o  where  i©’  is  the  intercept  at  ioo  =  0). 
Now,  having  k",  the  hydrogen  surface  coverage  0s  can  be  obtained  from  Equation  (6a), 
and  k3  can  be  obtained  from  Equation  (S). 


/ 


Thus,  Equation  (6c)  or  (6d)  quantitatively  illustrates  the  effect  of  U  on  the  h.e.r. 
The  effect  of  neglecting  ioo  is  to  eliminate  k”  in  equation  (6d)  so  that  one  only  obtains  the 
value  of  io'. 

Thus,  once  a  is  determined,  k",  ks  and  io’  can  be  determined  by  regression 
analysis  of  Equations  (5)  and  (6d),  if  the  plots  of  Uo  vs  (Equation  (5))  and  fi^^  vs  Uo 
(Equation  (6d))  are  linear.  An  absence  of  linearity  in  thd  experimentally  determined  plots  of 
Equations  (5)  and  (6d)  would  be  an  indication  that  at  least  one  of  the  three  assumptions 
mentioned  above  is  not  met,  e.g.,  steady  state  was  not  achieved  in  the  experiments,  the 
backward  reaction  of  the  discharge  step  is  significant  or  the  discharge-recombination 
process  is  not  operating.  Finally,  from  the  measured  values  of  ic  and  q,  0s  vs  q  can  be 
plotted  using  Equation  (6).  The  equilibrium  surface  coverage  of  hydrogen  (0e)  can  be 
determined  by  extrapolating  this  plot  to  q  =»  0.  Then,  from  Equation  (1) 


(7) 


a  ^e 
H 


-oe'V/rtI 


and 


(7a) 


io  =  io'  (1  -  0c)  (8) 

Comparison  of  the  discharge  rate  constant  and  the  recombination  rate  constant  can  provide 
some  insight  into  the  reaction  mechanism.  The  discharge  reaction  rate  is  potential 
dependent,  i.e.,  as  the  cathodic  overvoltage  (q)  is  increased,  the  potential  dependent 
reaction  rate  constant  (kirj)  will  be  increased;  whereas,  the  recombination  reaction  rate 
constant  depends  solely  on  physico-chemical  processes,  e.g.,  the  surface  diffusion  of  Had 
atoms.  Thus,  it  is  obvious  that  the  potential  range  has  to  be  considered  in  order  to  compare 
these  rate  constants  for  mechanistic  analysis.  From  Equation  (1)  the  potential  dependence 
of  kiTj  may  be  taken  to  be  of  the  form. 
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=  k.e-*'’’  (9) 

hi  1 

One  of  the  most  common  mechanisms  of  the  hydrogen  evolution  reaction  on  iron  is 
the  coupled  discharge-recombination  mechanism.  It  is  considered  that  an  overall  process  is 
coupled  if  the  energies  of  the  activated  states  of  more  than  one  step  are  within  1  to  2 
kcal/mole  with  respect  to  their  initial  statesO^)  since  the  rate,  vj,  is  proportional  to 

e-(AG. /RT)  being  the  activated  energy  of  reaction  i),  the  discharge  and  recombination 

>1*  # 

reactions  can  be  considered  to  be  coupled  if  viA'3  =  (e'^G^/RT)  /  (g-AG^/RT)  _ 

AG*)/RT  -  e±i5(Xy600  «  12  or  1/12.  Therefore,  if  vi  and  V3  are  within  one  order  of 
magnitude  of  each  other,  then, 

0.1  k3  <  kj^  <  10  kj  (10) 

and  the  hydrogen  discharge  and  recombination  reactions  will  be  coupled. 

Equation  (10)  can  be  utilized  to  estimate  the  potential  range  where  the  reaction  will 

ocevu:  by  a  coupled  discharge-recombination  mechanism.  From  Equations  (9)  and  (10),  the 
lower  limit  of  the  potential  range  (T|^)  is  given  by; 


tig  =  [ln(10  ki/k3)]/aa 


(11a) 


and  the  upper  limit  of  the  potential  range  (T|^  is  given  by: 
=  ln(lc/l0k,)]]  /aa 


(11b) 


Thus,  conditions  under  which  a  coupled  discharge  recombination  mechanism  operates  are 
that  k3  -  kiT^  and  i«i  vs  ViT  and  fiti  vs  Uo  are  linear.  These  plots  provide  a  convenient 
check  of  the  mechanism  since  if  either  plot  is  nonlinear  a  coupled  discharge-recombination 
mechanism  is  ruled  out. 

The  basic  model  developed  here  will  be  known  as  the  I-P-Z  model,  in  future 
references.  Although  this  model  was  developed  considering  acid  solutions,  similar 
relationships  will  be  valid  for  alkaline  solutions.  The  main  difference  will  be  that 


i 


in  Equation  (1)  will  have  to  be  replaced  by  aHio  (=1)  hence 

io'  =  Fkj  =  F  k  Also,  ki  =  k jajj+  e'O^'^^F/RT  jn  Equations 

(9)  through  (1 1)  will  be  replaced  simply  by  kj  =  k  je‘“^*^F/RT 

Another  impOTtant  modification  to  be  stressed  here  is  the  correction  to  be  applied 
when  6  is  larger  than  ~0. 1,  as  pointed  out  eaiiier.  In  the  above  derivation  of  the  basic 
model,  for  simplicity  sake,  Langmuir  conditions(Al3)  were  considered.  But  deviations  will 
occur  when  6$  is  significant  and  so  ic  and  ir  have  to  be  corrected  by  utilizing  the  Frumkin- 
Temkin  relationship^^)  for  the  hydrogen  coverage  function.  Then  Equations  (1)  and  (2) 
become 


(12) 


(13) 


Here  y  refers  to  the  rate  of  change  of  the  apparent  standard  free  energy  of 
adsorption  with  coverage(13).  Modifications  of  the  basic  model  given  above  using  Eqs. 
(12)  and  (13)  are  presented  elscwhete(8.M), 

The  question  as  to  whether  or  not  the  intermediate  reaction,  termed  the  absorption- 
adsorption  reaction,  is  actually  in  local  equilibrium  under  all  circumstances  needs  also  to  be 
critically  examined.  This  seems  especially  true  in  the  absence  of  any  direct  experimental 
evidence  on  the  question  and  in  light  of  the  fact  that  very  high  fugacities  can  be  involved 
under  cathodic  polarization  of  the  electrode(l^^^.  This  point  is  addressed  below  in  the 
Discussion. 


RESULTS  AND  DISCUSSION 
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I.  Data  Analysis  Using  the  Basic  I-P-Z  Model. 

The  basic  I-P-Z  model  developed  above  is  tested  for  various  literature  data(23.6,l8) 
of  hydrogen  permeation  experiments  on  iron  of  different  purity  and  zone-refined  (99.995 
wt  %)  nickel.  Rgures  1  and  2  show  the  plots  of  log(ic)  vs  ^  and  log  (ioo)  vs  T|  for  these 
cases.  It  is  clear  that  the  plots  are  linear  and,  hence,  the  quadratic  equation  in  a  (Equation 
(A.5)  in  the  Appendix)  is  easily  solved  for  each  of  the  cases,  and  these  a  values  are  given 
in  Table  I. 

Figure  3  shows  the  Ue  vs  Vlr  plots,  which  are  all  linear  as  confirmed  by  various 
statistical  tests.  Three  of  the  plots  pass  through  or  close  to  the  origin  and  one  has  a 
relatively  large  non-zero  (positive)  intercept  By  performing  various  statistical  analyses, 
including  the  Student's  t  Test^^^l,  it  has  been  confirmed  that  the  non-zero  (positive) 
intercept,  for  the  data  of  Bockris  et  alS^,  is  significant  at  95%  or  even  at  99%  confidence 
levels.  The  basic  I-P-Z  model  is  strictly  valid,  however,  only  for  plots  that  pass  through 
the  origin.  Therefore,  for  the  three  cases  in  Hgure  3  where  the  Ue  vs  Vv  plots  are  linear  and 
pass  through  the  origin,  the  basic  I-P-Z  model  is  applicable  and  so  ic^  versus 
(Equation  (6d))  can  be  plotted.  Figure  4  shows  these  plots,  which  are  seen  to  be  linear, 
confirming  the  assumption  of  steady  state  [Le.,  Assumption  (4)].  But  due  to  the 

I 

extremely  small  values  of  the  slopes  in  Bg.  4  are  quite  small  making  the  lines  seem 
almost  horizontal.  F<v  the  plot  in  Fig.  3  widi  the  non-zero  intercept,  modificaticm  of  the 
basic  I-P-Z  caodel  is  necessary  and  this  is  presented  in  Section  n. 

Since  the  plots  in  Figs.  3  and  4  are  linear,  all  of  the  constants  k3,  io',  k",  k°  and  9e 
are  determined  for  the  data  in  References  3, 6  and  18  from  the  regression  analysis  of 
Equations  (5)  and  (6d).  A  value  of  Di  =  5  x  10*5  cm2  used  for  iron  in  these 
analyses  although  measurements  of  the  true  lattice  hydrogen  diffusivity  in  iron  are  not 
easily  made  as  shown  by  the  great  variadon  in  measured  values(29).  For  nickel,  Di  =  4  x 
10-10  cm2  s-l(l8)  was  used.  All  of  these  values  are  tabulated  in  Table  I. 
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The  accuracy  levels  of  the  computed  values  of  the  transport  and  kinetic  constants 
and  coefficients  in  Table  I  are  primarily  limited  by  the  accuracy  of  the  diffusivity  and  the  a 
values.  These  results  were  confirmed  by  various  statistical  analyses  (including  the  F  test 
and  the  Student's  t  test)  in  that  the  confidence  levels  of  the  estimated  parameters,  (i.e., 
slopes  and  intercepts),  from  the  regression  analyses  of  Equations  (5)  and  (6d),  were  better 
than  that  possible  for  the  diffusivity  and  a  values.  Even  a  slight  error  in  a  can  cause  a 
large  scaner  in  the  ice®“n  vs  "n  plot,  since  a  appears  in  the  exponential  term. 

Since  k"  is  known,  0$  values  are  obtained  by  using  Equation  (6a)  for  various  iee  (or 
equivalently  Ti)  values.  0s  vs  r\  plots  are  given  in  Figure  5,  which  shows  that  the  hydrogen 
surface  coverages  are  quite  small  in  all  of  the  cases,  i.e.,  0s  <  0.  IS.  Although  only  a  few 
experimental  determinations  of  0$  exist  in  the  literature,  values  for  comparison  with  those 
in  this  paper  for  the  0. 1  N  H2SO4  solutions  have  yielded  similar  0s  values,  i.e.,  in  the 
range  of  0.01  to  0.1(10,21).  More  (direct)  experimental  measurements  of  0s,  however,  are 
needed  for  checking  the  model. 

As  shown  above,  this  model  is  able  to  determine  the  transfer  coefficient  oc,  as  per 
Equation  (A.S)  of  the  Appendix  and  it  is  a  function  of  both  the  inverse  Tafel  slope,  Sc,  and 
the  slope  ($«),  dln(i«,)/dTi.  Table  I  shows  that  the  a  values  fall  on  either  side  of  0.5  for  the 
acid  solutions  with  the  data  of  Kato  et  al.O^  being  tire  furthest  firom  0.5.  It  is  to  be  noted 
that  the  different  iron  data  were  for  sanqiles  having  varying  impurities,  membrane 
thicknesses  and  pretreatments.  The  tenttency  for  surface  segregation  (e.g.,  during  prior 
heat  treatment)  is  greater  the  less  ptire  is  the  iron,  and  the  tendency  fra*  impurity  adsorptirat 
(fiom  the  electrolyte  during  the  permeatirai  expoiment)  is  greater  the  thicker  is  the  sanqile 
because  of  the  longer  permeation  time.  Kato  et  used  99.98  wt  %  iron  (0.004%  C, 
0.006%  Si  and  0.002%  P,  S  and  Al)  membranes,  2mm  thick  and  0.1  N  H2SO4  charging 
solution,  and  the  charging  surface  was  mechanically  polished  prior  to  the  hydrogen 
charging  experiment.  Bockris  et  used  Armco  iron  (99.9  wt  %)  membranes,  0.77  mm 
thick,  0.1  N  H2SO4  and  the  charging  surface  was  etched  as  tiie  final  sample  preparation 
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stq).  Zamanzadeh  et  used  Ferrovac  E-iron  (the  most  pure  of  the  three)  membranes, 
0.365  mm  thick,  alkaline  solution  (O.IN  NaOH)  and  the  membranes  were  annealed  at 
1273  K  in  an  argon  back-filled  evacuated  apsw  \  In  addition,  the  surface  pH  values  were 
probably  different  for  the  three  sets  of  data  even  for  the  O.IN  H2SO4  soluticms  in  the 
absence  of  a  buffer.  Considering  all  of  these  aspects,  it  is  not  surprising  to  see  the 
variability  of  the  computed  values  for  the  three  sets  of  iron  data  in  Table  I.  From  Equaticm 
(1),  an  explanation  of  why  the  io  values  for  iron  are  different  in  each  case  would  be  based 
on  pH  and  factors  affecting  both  k^  (such  as  impurity  adsorption)  and  a. 

Besides  surface  conditions  such  as  the  presence  of  oxide  films,  other  factors,  such  as 
impurity  adsorption,  surface  segregation  and  surface  enrichment  that  is  caused  by  selective 
dissolution,  can  change  the  kinetics  of  the  surface  reactions,  including  the  absorption- 
adsorpdon  step  and  hydrogen  desorption.  The  effect  of  one  or  more  of  these  facttxs  could 
account  for  the  variation  of  k"  (besides  the  very  important  thickness  effects,  as  described 
later)  and  kj  (recombinaticm  constant)  values  for  the  iron  data  in  Table  1 . 

As  pointed  out  in  the  last  section,  the  potential  dependent  discharge  reactimi  rate 
coefficient  (kit\)  depends  on  the  overvoltage  (r|).  At  very  low  Tl  values,  the  discharge 
reaction  will  be  slower  than  the  recombinatitm  reaction  and  at  moderate  to  high  overvoltages, 
they  will  be  coupled.  In  general,  for  iron  (refer  to  Table  I)  hydrogen  charging  at  an 
oveipotential  between  -400  and  -900  mV  will  involve  a  coupled  discharge-recombinadon 
mechanism.  For  nickel,  this  overpotendal  range  is  between -440  and -650  mV.  However, 
the  overpotendals  employed  in  the  experiments  were  smaller  than  these  ranges  (see  Figures 
1, 2,  or  5),  i.e.,  the  overpotendals  correspond  to  the  discharge  reacdon  being  slower  than  or 
about  as  fast  as  the  recombination  reaction,  consistent  with  the  low  values  of  0$  in  all  of  the 
cases  (Hg.  5).  These  observations  may  provide  a  clue  in  understanding  the  mechanism  in 
the  fanuliar  industrial  problem  of  the  enhancement  of  hydrogen  permeatimi  in  the  presence  of 
poisons  like  H2S,  in  that  a  reduction  in  k3  (e.g.,  a  reduction  in  surface  diffusivity  of  Had)  or 
increase  in  k^  or  both  can  push  the  reaction  towards  a  fully  coupled  mechanisia  In  this  way. 


due  to  the  presence  of  poisons,  such  as  H2S,  the  hydrogen  coverage  is  increased  and  since 
the  hydrogen  diffiisivity  is  relatively  high  in  iron,  the  hydrogen  permeation  is  greatly 
enhanced.  A  more  advanced  study  of  these  aspects  in  the  case  of  enhanced  hydrogen  entry 
due  to  poisons  is  shown  in  the  modified  I-P-Z  model(8.l4)  taking  into  account  the  Frumkin- 
Temkin  corrections(*3)  for  hydrogen  coverage, 
n.  A  Modified  I-P-2  Model  for  Selvedge  Reactions. 

The  fact  that  a  statistically  si^iificant  non-zero  (positive)  intercept  was  obtained  in 
Figure  3,  for  the  data  of  Bockris  et  rz/.(2),  poses  an  interesting  question:  Is  Assumption  (2) 
valid  in  all  cases,  i.e.,  is  the  adsorption-absorption  reaction  in  local  equilibrium?  As  pointed 
out  earlier,  in  the  absence  of  any  direct  expoimental  evidence  on  the  question  and  in  the  light 
of  the  fact  that  very  high  fugacities  can  be  involved  under  cathodic  polarizatitxi  of  the  metal 
electrode(’^^'l‘^  there  is  good  reason  to  consider  other  possibilities  for  the  intermediate 
reaction  that  follows  the  hydrogen  discharge  reaction.  In  terms  of  quantitatively  treating 
hydrogen  permeation  data  that  contain  positive  or  negative  intercepts,  it  means  considering 
other  intermediate  reactions  that  may  themselves  produce  concentration  profiles  in  the  near 
surface  (V  selvedge  region  that  differ  the  concentration  profile  in  the  bulk  metal 
(established  by  the  bulk  diffusion  of  hydrogen).  Three  different  possible  cases  of  H 
concentration  profiles  are  sketched  in  Figure  6. 

Figure  6(a)  depicts  the  assumed  condition  for  the  basic  I-P-Z  noodel,  i.e.,  the 
conditiem  of  local  equilibrium  of  the  adscaption-absorptiem  reaction.  In  this  case  U,  is 
proportional  to  VW « tuid  the  intercept  (Cg)  of  tiie  Uo  vs  Vir  plot  is  zero.  Furthermore,  fev 
equilibrium  of  the  adsorption-absorption  reaction,  the  physical  situation  is  that  the  absorbed 
state  is  located  in  the  very  near  vicinity  of  the  adsorbed  state,  i.e.,  implicitly  an  atomic 
dimension  below  the  surface. 


The  other  two  cases  in  Figure  6  (i.e.,  6(b)  and  6(c))  consider  a  different  physical 
situation  that  is  referred  to  as  a  selvedge  in  the  gas-solid  literature(^).  The  selvedge  in  the 
context  of  electrochemical  H  charging  can  be  defuied  as  an  intermediate  reaction  layer,  more 
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than  an  atomic  layer  thick,  having  continuity  with  the  bulk  substrate.  Essentially,  the  concept 
of  a  selvedge  in  this  paper  will  be  used  to  describe  the  specific  subsurface  region  of  hydrogen 
interaction  with  the  metal.  The  situation  in  Hgure  6(c)  differs  fix>m  that  in  Figure  6(a)  only 
in  that  the  hydrogen  diffusivity  in  the  selvedge  is  greater  than  that  in  the  bulk  substrate  (if  the 
hydrogen  diffusivity  were  lower  in  the  selvedge  than  in  the  bulk,  the  slope  would  be  of  the 
same  sign  but  steeper  than  that  of  the  bulk).  The  situation  in  Figure  6(c)  might  depict,  for 
example,  a  surface  segregated  condidon.  In  this  case  the  intercept  (Cg)  of  the  i«o  vs  VIr  plot 
will  be  negative. 

The  question  still  remains  of  the  case  where  the  intercept  (Cg)  of  the  ioo  vs  plot  is 

positive  as  in  the  data  of  Bockris  et  al.  (Figure  3).  For  a  pure  metal,  this  can  only  be 
achieved  by  a  diffusionless  process  such  as  a  situation  where  the  hydrogen  atom  actually 
penetrates  into  the  substrate  a  finite  distance,  li,  (of  atomic  dimension)  immediately  after  the 
proton  discharge  step.  This  possibility  has  not  been  ruled  out  on  theoretical  grounds  and  the 
above  mentioned  positive  intercept  in  the  data  in  Reference  2  indicates  that  some  type  of 
difhisionless  process  may  be  involved.  If  hydrogen  does  penetrate  into  the  lattice  driven  by 
the  energy  associated  with  the  proton  dischar^  step,  the  initial  hydrogen  subsurface 
concentration  profile  should  be  similar  to  that  obtained  during  ion  implantation,  in  that  a 
maximum  in  concentration  occurs  at  some  distance  beneath  the  substrate  surface  where  the 
penetrating  atoms  come  to  rest  The  normal  hydrogen  diffusional  process  would  ccMiunence 
immediately  following  this  nondiffusional  step  with  most  of  the  hydrogen  atoms  diffusing 
back  towards  the  input  surface  since  this  surface  is  much  closer  to  the  concentration 
maximum  than  is  the  other  (exit)  surface  of  the  membrane.  Since,  during  ongoing  hydrogen 
discharge,  both  processes  (penetration  and  diffusion)  may  occur  simultaneously,  the  profile 
for  steady  state  hydrogen  permeation  could  be  as  shown  in  Figure  6(b).  Thus,  the  intercept 
(cg)  in  the  Uo  vs  Vir  plot  (Figure  3)  is  positive. 


In  order  to  take  care  of  the  various  possible  cases  of  concentration  profiles  (Figure 
6(b),  where  Cg  >  O  and  Figure  6(c),  where  Cg  <  O)  the  basic  I-P-Z  model  has  to  be  modified 
slightly,  as  follows: 

Cg=  Ci  -  Cs  (14) 

where  Ci  =  the  subsurface  hydrogen  concentration  at  the  selvedge  boundary  and  Cg  is  the 
concentration  difference  in  the  selvedge,  which  is  assumed  to  be  independent  of  ic',  ie., 

Cs  =Ci  -Cg  (14a) 

Equation  (3)  will  be  modified  as: 

i«zF(Di/L)ci  (14b) 

From  Equations  (14a)  and  (14b), 


c. 


Km' 


(15) 


wheiob  =  L/(FDi). 

From  Equations  (2),  (4)  and  (15), 

i.  .  +  ('/=) 


(16) 


Similarly,  using  Equations  (1),  (4)  and  (15),  Equation  (6d)  will  be  noodified  (by 
considering  the  selvedge)  as: 


»c  e 


aocq  __ 


bi. 


k"; 


(17) 


Utilizing  Equations  (16)  and  (17),  the  data  of  Bockris  et  alP-)  can  be  analyzed  as  was 
done  above  for  the  other  three  cases.  The  modified  plot,  ic  e^^n  vs  fioo  -  (cg/b)]  is  shown  in 
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Hgure  7.  The  relevant  quantities  utilized  for  this  analysis  and  the  resulting  computed 
parameters  are  shown  in  Table  I,  alongside  the  results  for  the  other  three  cases.  It  is  to  be 
noted  in  Table  I  that  the  value  of  Cg  is  quite  small  in  this  case  and  so,  neglecting  the  value  of 
Cg  in  calculations  of  the  various  transport  and  kinetic  parameters  would  not  introduce  serious 
errors. 

A  modified  version  of  the  reaction  sequence  used  in  the  derivation  of  the  basic  I-P-Z 
model  is  shown  in  Figure  8,  where  the  selvedge  reaction  is  also  included, 
m.  The  Limiting  Conditions  for  Diffusion  Control  Versus  Interface  Control. 

From  Equation  (4), 


Substituting  for  Cs  (=  p^i-)  fro®  Equatitm  (3)  into  Equation  (18a)  and  simplifying: 


tL=(  +  ( _ ! — ) 

FOikabsOs  Fkabs8s 


(18b) 


Equation  (1 8b)  describes  the  same  relationship  as  derived  by  Kim  and  Wilde^O)  and  by 
Ateya  et  alS^  i).  As  Ateya  et  a/.O  i)  have  pointed  out,  there  are  two  iroptntant  limiting  cases  of 
Equation  (18b)  or  equivalently  Equation  (4). 

Case  I.  For  thick  membranes,  or  low  values  of  Di,  i.e.  under  conditions  of  pure  diffusion 
control,  kads  and  hence  Equation  (4a)  will  reduce  to: 


As  pointed  out  before,  k'  is  usually  termed  as  the  equilibrium  abscnption-adsorption 
constant. 


It  can  then  be  easily  shown  that  Equation  (18b)  reduces  to: 
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1  _  kads  _  1 

^  ~  FkabsOs  "  Fk’0s 

Rearranging,  one  has 


(19b) 


u.  =  F^k’es  (19c) 

Equation  (19c)  is  a  familiar  relationship.  In  other  words,  the  assumption  of  the  intermediate 
reaction  rate  constant  being  thickness  independent  can  be  valid  for  thick  membranes  or  low 
values  of  Di. 

Case  TI.  For  very  thin  membranes,  or  high  values  of  Di,  or  low  values  of  kabs» 
and  hence  Equation  (4a)  will  reduce  to: 


It  can  then  be  easily  shown  that  Equation  (18b)  reduces  to: 

i.  ^ 

FkabsQs 

io*  =  Fkabs^s 


(20a) 


(20b) 


(20c) 


Equation  (20c)  is  the  same  relationship  as  that  of  Ateya  etalS^^\  As  pointed  out  by  them,  i«, 
becomes  thickness  independent  and  the  hydrogen  entry  will  be  interface  controlled. 

IV.  Evaluation  of  the  Surface  and  Subsurface  Kinetic  Properties  of  a  Metal  Electrode  Using 
the  I-P-Z  Model. 

The  H  adsorption  and  absorption  rate  constants  describe  the  surface  and  subsurface 
kinetics  for  the  metal  electrode.  These  kinetic  properties  could  eventually  be  utilized  to 
characterize  the  surface  and  subsurface  states  of  a  metal. 
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From  Equation  (4a),  the  following  important  general  relationship  can  be  obtained: 


1 


(21) 


The  hallmark  of  Equation  (21)  is  that  by  plotting  vs.  p  one  can  obtain  some  of 
the  ever-sought  surface  properties,  such  as  kabs  and  k«Hc.  If  the  ^  vs.  ^  plot  is  linear,  then 
Equation  (21)  is  applicable,  which  implies  that  the  adsorbed  hydrogen  is  the  same 
intermediate  which  takes  part  in  both  the  evoludcm  and  absorption  steps,  as  has  been 
proposed  by  Bockris  and  others(2.4,23,24)  oygj.  jjjg  (}ecades  and  generally  accepted. 

Figure  9  shows  the  plot  of  vs.  ^  for  iron  data  (as  per  Table  I).  It  shows  the 
correct  trend,  as  expected  from  Equation  (21),  illustrating  the  validity  of  k"  values  obtained 
by  the  application  of  the  I-P-Z  oxxlel  to  the  permeation  data  of  various  iron  membranes.  But 
it  is  to  be  noted  that  vs.  ^  is  not  linear,  even  after  considering  the  erTW  bands.  This, 
however,  does  not  imply  that  the  adsorbed  and  absorbed  hydrogen  are  not  the  same 
intermediates,  since  in^urity  levels  of  die  iron  membranes  as  well  as  die  electrolytes  are 
different  for  the  three  sets  of  iron  data  and  hence  the  kabs  (and  even  DO  and  kads  values  would 
be  different 

V.  Advantages  of  the  I-P-Z  Model 


The  I-P-Z  model  has  the  following  advantages. 

(1)  It  enables  calculation  of  surface  coverages,  reaction  rate  constants  and  the 
potential  range  in  which  a  particular  reaction  mechanism  will  occur,  by  utilizing  w  for  the 
determination  of  k",  etc.  as  noted  earlier  in  conjunction  with  Equations  (6a)  through  (6d). 

(2)  It  is  able  to  explain  physically  and  mathematically  all  the  different  processes 
involved  in  a  self  consistent  manner.  Also,  it  is  quite  flexible  to  nuxlifications  that  will  take 
care  of  changing  physical  and  chemical  conditions.  That  is,  the  basic  I-P-Z  model  will  serve 
as  a  building  block  for  more  sophisticated  models. 
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(3)  The  I-P-Z  model  in  conjunction  with  the  derivations  of  Ateya  et  for  the 
first  time  focuses  on  the  importance  of  the  membrane  thickness  and  the  limiting  conditions 
for  diffusion  control  versus  interface  control. 

(4)  It  has  a  unique  capacity  to  quantitatively  evaluate  the  ever-sought  surface  and 
substuface  kinetic  properties,  i.e.,  the  H  adsorption  and  absorption  rate  constants. 

(5)  It  delineates  the  roles  of  electrochemical  kinetics  (discharge  reaction  -  proton 
tunneling),  possible  unique  surface  layer  kinetics  (selvedge  reacticm)  and  physico-chemical 
kinetics  (permeation  and  recombination  reactions).  Strictly  speaking,  all  of  the  rate  constants 
(except  that  for  the  discharge  reaction)  should  be  potential  independent  and  this  model 
recognizes  that  From  these  simple  concepts,  the  model  advances  that  the  potential  range  in 
which  a  certain  mechanism  will  be  operating  is  bounded.  This  is  true  for  the  coupled 
discharge  -  recombination  mechanism,  analyzed  by  the  model.  Comparison  of  the  various 
calculated  rate  constants  help  to  verify  the  mechanism  and  the  validity  of  the  assumptions. 
These  ate  further  checked  by  the  linearity  oi  the  two  plots,  u.  vs  Vir  and  fit]  vs  u’.  In  short, 
the  iixxlel  can  serve  as  a  diagnostic  scheme  to  analyze  consistency  between  the  observed 
relationship  and  the  proposed  mechanism  by  calculating  the  rate  constants  and  the  critical 
potential  ranges  from  the  data  base. 

(6)  The  model,  when  adequately  modified,  has  the  capacity  to  explain  and  predict  the 
behavior  of  poisons,  such  as  H2S,  enhancing  permeation  but  in  many  cases  decreasing  the 
hydrogen  overvoltage  at  a  constant  charging  density.  This  will  be  considered  in  detail  in  an 
upcoming  papei<^),  where  the  model  is  irxxlified  by  incorporating  the  Friimkin-Temkin 
correction. 

VI.  Limitations  of  the  I-P-Z  Model 

(1)  The  model  cannot  be  used  in  its  present  form  when  the  backward  reactions 
become  important,  i.e.,  at  small  hydrogen  oveipotentials  (T)  <  RT/F) 

(2)  The  possibility  of  the  occurrence  of  the  selvedge  reaction  and,  hence,  the  value  of 
li  have  to  be  experimentally  and  theoretically  investigated. 
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(3)  The  actual  potential  dependence  of  kir^  (Equation  (9))  as  well  as  the  proper 
criteria  for  finding  the  potential  range  in  which  a  mechanism,  such  as  the  coupled  discharge- 
recombination  reaction  operates,  have  to  be  evaluated. 

Notwithstanding  the  above  limitations  of  this  model,  it  can  provide  some  important 
information  about  surface  coverages,  rate  constants  and  opwating  mechanism.  The  concept 
of  the  selvedge,  if  proven  by  specially  designed  experiments,  has  the  capacity  to  provide  a 
detailed  physical  picture  for  explaining  many  of  the  complex  phenomena  occurring  during 
hydrogen  entry  and  embrittlement  in  many  materials. 

CONCLUSIONS 

1.  A  model  has  been  developed  which  for  the  first  time  takes  into  account  the  effect  of  the 
hydrogen  permeation  flux,  U#,  on  the  hydrogen  evolution  reaction  during  electrochemical 
hydrogen  charging  of  materials.  This  model  gives  the  critically  important  values  of  the 
absorption  (kabs)  and  adsorption  (kads)  rate  constants,  discharge  and  recombination  rate 
constants  and  hydrogen  coverages  from  hydrogen  permeation  data.  In  addition,  an  analytical 
method  in  the  model  facilitates  the  determination  of  the  operating  mechanism  of  the  hydrogen 
evolution  reaction. 

2.  The  model  also  incorporates  the  thickness  and  diffusivity  factors  of  the  membrane  into  the 
absorption  and  adsorption  rate  constants  in  order  to  address  and  compare  the  various 
electrochemical  and  physico-chemical  reactions. 

3.  A  few  of  the  literature  data  of  hydrogen  permeation  in  iron  and  nickel  are  successfully 
analyzed  with  this  model. 

4.  Further  experiments  and  analyses  need  to  be  performed  in  order  to  check  the  validity  of 
the  assumptions,  and  the  role,  if  any,  of  the  selvedge  reaction. 
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UST  OF  SYMBOLS 


a  a  constant,  F/RT,  (volts)' ^ 

af{+  hydrogen  ion  activity,  dimensionless 

aH20  activity  of  water  in  the  electrolyte,  dimensionless 

A  amperes 

b  a  constant,  L/FDh,  (A.L)‘^ 

Cg  concentration  difference  in  the  selvedge,  Ci-Cg,  moll^^ 

CH+  hydrogen  ion  concentration,  lO'P*^,  nnol 

Ci  subsurface  hydrogen  concentration,  nwl 

Cs  surface  hydrogen  concentration,  mol 

Di  hydrogen  diffusion  coefficient  in  the  selvedge  (transition  layer), 
Di  lattice  hydrogen  diffusion  coefficient,  S'* 

1^  equilibrium  potential  for  the  h.e.r.,  mV  vs  SHE 

F  Faraday  constant,  965(X)  C  (g-eq)** 

fin  a  variable,  ice*®*!,  A 

AGi*  activation  energy  of  reaction  i,  J(g-mol)'* 

ic  charging  flux  or  current  density,  A 1^'^ 

ioo  steady  state  permeation  flux,  A 

ioo'  a  variable,  w-Cg/b,  A 

io  exchange  current  density,  A 

io’  io/(l-ee).AL-2 

if  steady  state  evolution  flux,  A  L*^ 

i*abs  absorption  rate  constant,  mol  (L^s)'* 

kads  adsorption  rate  constant,  L  s'* 

k '  equilibrium  absorption-adsorption  constant,  mol 

k ' '  thickness  dependent  absorption-adsorption  constant,  mol 

k  1 0  discharge  reaction  rate  constant,  mol  (L^s)'* 


ki  discharge  reaction  rate  coefficient,  nK>l  (L^s)*^ 
klTi  potential  dependent  discharge  reaction  rate  coefficient,  mol  (L^s)* 
k2  selvedge  (intermediate)  reaction  rate  constant,  L.s*^ 

k-2  reverse  (intermediate)  reaction  rate  constant,  L.s'  * 
k3  recombination  reaction  rate  constant,  nool  (L^)'  ^ 

L  membrane  thickness,  L 

L  unit  length 

li  selvedge  depth,  L 

R  gas  constant,  8.314  J  (g-mol  K)*l 

Sc  slope  of  In(ic)  vs  t],  decade  V** 

Set  slope  of  7]  vs  log  ic,  V  decade'^ 

Soe  slope  of  ln(i«o)  vs  ti,  decade  V** 

Soot  slope  of  ti  vs  log  ioo,  V  decade** 

T  temperature,  K 
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Greek  Symbols 

a  transfer  coefficient,  dimensionless 

Tj  overvoltage,  V 

Tic  cathodic  overvoltage,  V 

Tjg  lower  limit  of  Tic,  V 

upper  limit  of  T|c,  V 

6s  surface  coverage,  dimensionless  • 

de  equilibrium  surface  coverage,  dimensionless 

Y  gradient  of  the  apparent  standard  free  energy  of  adsorption  with  hydrogen  coverage, 
J  (g-mol)-l 
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APPENDIX 

Determination  of  the  Transfer  CoeffidenL  a; 


In  previous  determinations  of  oc,  6$  was  assumed  to  be  close  to  zero  and  neglected 

dTi 

which  resulted  in  the  familiar  relationship  between  a  and  the  Tafel  slc^  )= 


2  3  RT 

- :  However,  in  keeping  with  the  thrust  of  this  papa*,  Le.,  taking  into  account  the 

aF 

important  effect  of  6$  which  is  potential  dependent,  a  needs  to  be  rigorously  derived  as 
follows.  Equation  (6c),  in  the  text,  can  be  rewritten  as: 


ioo  =  - 


+  (k‘/b) 


(A.1) 


Differentiating  Equation  (A.1)  with  respect  to  the  oveipotential,  Ti,  (assuming  a  and  k" 
are  independent  of  ‘r\)  one  obtains 

(di JdTi)  =  -( kXbi'o)  )  [aodce*®^  +  c*®^  (dic/dn)]  (A.2) 

and 

[d  ioo/dTl^)  =  -  (k  /(^^o))e*®^^d^iydTi^]  +  2aa(di</dT|)  +  (aa]^]  (A.3) 

Dividing  Equation  (A.3)  by  Equation  (A.2),  and  simplifying,  one  obtains 

+  (o/a)  ^-|(d^i«/dTi^]/(di«/d‘n]|  +  {(2/ic)(dio/dTi)}j  + 

(l/a^)  [(l/ic)(dVdn^)  -  (l/ic)(dic/dTi){(d^i  Jdn^)/(diJdn)}]  =  0  (A.4) 

From  plots  of  ic  vs  ti  and  Uo  vs  ti,  the  above  differentials  (in  Equation  (A.4))  can  be 
determined  and  then  the  quadratic  equation  in  a  can  be  solved.  If  the  measured  In  (ic)  vs  T| 
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and  In  (Uo)  vs  t|  plots  are  linear,  this  con&ms  that  a  is  independent  of  t],  and  Equation 
(A.4)  can  be  simplified  as: 


a2  +  ((2sc-sj/a)a  +  (sc(Sc-sJ/a^)  =0 


(A.5) 


where  Sc  =  dln(ic)/dTi  and  s«,  =  dln(ieo)/dTi.  The  prq>cr  value  of  a  is  chosen  from  the 
soludon  of  Equation  (A.5)  so  as  to  obtain  non-negative  values  of  Jo'  and  k"  in  Equation 
(A.1). 
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TABLE  I 

Analysis  of  Literature  Data  (References  as  shown^ 


Value 


Parameter 

Ref[3] 

Refl6] 

Reftl8] 

Ref[2] 

Material 

PureFe 

Ferrovac  Fe  I 

Zone-reEned  Ni  Armco  Fe 

Di  (cm^  s'^) 

5  X  10*5 

5x10-5 

4x10-10 

5x10-5 

L(mm) 

2 

0.37 

0.08 

0.77 

b  (mole(A-  cm*!)) 

0.04 

0.01 

200 

0.02 

pH 

1 

13 

1.3 

1 

a 

0.38 

0.45 

0.56 

0.56 

io  =  io’(Acm-2) 

3  X  10*6 

1x10-6 

4x10-6 

6x10-7 

kj  (cm  s*^) 

lx  10-7 

0.2 

2x10-7 

2  X  10-8 

ki  (mole  (cm^  s)-!) 

4x10-11 

1  X  10-11 

4x10-11 

6x10-12 

k3  (mole  (cm^-s)*^) 

3x10-6 

1x10-7 

5x10-6 

7  X  10-11 

k"  (molecm*3) 

2x10-6 

2x10-7 

8x10-3 

1x10-6 

Cg  (mole  cm*3) 

4x10-8 

0c 

0,004 

0.01 

0.003 

0.003 

Tli(mV) 

-610 

-380 

-440 

-430 

Tic  (mV) 


-920 


-640 


-650 


-640 
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HGURE  CAPTIONS 

1 .  Various  current-potential  data  for  the  hydrogen  evolution  reaction  on  iron  and  nickel 
electrodes. 

2.  Various  hydrogen  permeation  data  for  hydrogen  charging  of  the  iron  and  nickel 
electrodes  in  Fig.  1. 

3.  Plots  of  the  hydrogen  permeation  current  versus  square  root  of  the  recombination 
current. 

4.  Plots  of  the  hydrogen  charging  function  (icC®*^)  versus  the  steady  state  hydrogen 
permeation  (ioo),  for  iron  at  pH  1,  iron  at  pH  13  and  nickel  at  pH  1.3. 

5.  Plots  of  the  fractional  hydrogen  coverage  (8$)  as  a  functimi  of  the  overpotential  (t|) 
obtained  by  applicatitxi  of  the  model  to  the  data,  (a)  q  vs  8$.  (b)  r\  vs  log  8$. 

6.  Sketches  of  the  three  different  cases  of  H  concentration  profiles,  expected  for  a  general 
(physicsd)  intermediate  reaction. 

7.  Plot  of  the  hydrogen  charging  function  (ice*“n)  versus  the  hydrogen  permeation 
function  (iao  -  Cg/b),  for  iron  at  pHl. 

8.  Schematic  showing  hydrogen  discharge,  recombination,  permeation  and  selvedge 
reactions,  specificdly  for  the  case  of  Hg.  6(b). 

9.  Plot  of  the  absorption-adsoiption  constant  (k**)  as  a  function  of  the  membrane  thickness 
(L)  for  various  data  on  iron. 
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